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ABSTRACT 

Historically, studies on decorrelation scales have been conducted in the deep 
ocean waters. As the Navy shifts its interest toward the less understood shallow 
water regions, decorrelation scales need to be computed in order to use formerly 
deep water models such as the Optimum Thermal Interpolation System (OTIS) for 
shallow water regions such as the Yellow Sea. 

A data set containing over 35,000 temperature profiles from 1929 to 1991 
was obtained from the Naval Oceanographic Office’s MOODS data set. The winter 
and summer seasons provide realistic results. Winter has the smallest decorrelation 
scales of all the seasons, approximately 15 days and 165 km. Summer shows that 
there are different decorrelation scales between the surface and at depth. The surface 
has scales of 12.3 days and 251 km while at depth the scales are approximately 16.5 
days and 163 km. An observational sampling network design is suggested for future 
sampling of the region. Spring and fall provide mixed results which may be due to 
the irregularities in time and space of the data set or to the very complex forcing 
mechanisms found in the region. Overall, this study gives a ground work for better 
refinement of decorrelation scales and thus, the ability to assess the conversion of 
deep water models to shallow water regions. 
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I. INTRODUCTION 



Various ocean systems such as fronts, eddies and water masses have different 
temporal and spatial correlation scales (Figure 1.1). These scales feature a system’s 
life span and spatial extent both horizontally and vertically. Horizontal scales are 
often divided into north-south and east-west length scales and are represented by a 
variety of parameters. For example, Provost and Le Traon (1993) attempted to 
describe the Brazil-Malvinas Current confluence confines through temporal and 
spatial scales using variable altimetric measurements. 

Horizontal scales for temperature fields in the deep ocean basins have been 
extensively studied. Each ocean has its own general scales with variations in specific 
regions. White et al. (1982) identified spatial correlation scales in the western 
Pacific of about 600 km in the tropics (south of 17.5°N) and 300 km in the subtropics 
(north of 17.5°N) while Ozsoy et al. (1989) found the spatial scale to be 200-250 km 
in the Levantine Basin of the Mediterranean Sea. Historically, most of the studies on 
correlation scales have been conducted in deep ocean waters; few have been 
conducted on the continental shelf. There are several reasons for this occurrence 
with the main one being that deep water physics are better understood. Another, yet 
not an obvious reason, is that in past years the US Navy has put its emphasis on deep 
water operations. Therefore, most of the Navy’s research and modeling assets have 
been directed toward deep water applications. This poses many challenging 
problems as the Navy shifts its interests toward the less understood shallow waters 
regions. Naval models such as the Optimum Thermal Interpolation System (OTIS) 
that have been designed for deep water usage need to be modified for shallow water 
use. 

Deep ocean temperature profiles tend to follow the pattern of a mixed layer, 
a thermocline and a deeper layer slowly decreasing in temperature with depth. 
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Figure 1.1 Approximate Ranges of Horizontal and Temporal Scales of 
Selected Ocean Variabilities (fromMooers and Robinson, 1981). 
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Shallow water does not consistently mimi c this pattern. It may range from the 
classical profile or be completely isothermal. The shallow environment influences 
the temperature distribution with some parameters like bathymetry, river run-off and 
tides, often omitted in deep water, becoming increasingly important in shallow 
water regimes. Spatial scales are strongly controlled by bathymetry. In shallow 
water, besides scales associated with seasonal changes, some other time scales also 
become important. The forcing of internal waves (minutes), tides (hours-day), 
meteorological forcing (hours-weeks), freshwater runoff (days-months), and 
interactions with slope or deep water (days-weeks) can have a great impact on 
temporal variability (Haeger, 1991). 

The Yellow Sea was selected for temperature horizontal correlation analysis 
based on two main considerations. The first is that the Yellow Sea demonstrates 
shallow water characteristics for which a data set of 50,000 profiles exists. The 
second consideration is from a military standpoint. The region is a potential “hot 
spot” and needs special study. Modification of the OTIS model for the Yellow Sea 
can be used for other marginal seas. To do so, there must be an accurate climatology. 
Since this is a new region of emphasis, as stated earlier, a climatology has not been 
built from the Naval resources. Thus, there are two tasks in this study, (1) 
constructing a seasonal climatology from the data base and (2) computing horizontal 
and temporal scales and their vertical variation. This work will provide useful 
information for the OTIS model, since it needs temporal and spatial scales prior to 
running the model. After the climatology has been established, temperature 
anomalies are computed by subtracting the climatology from the observations. The 
anomalies are then paired and placed into bins depending on their separation in time 
and space. A correlation coefficient is then calculated from all the anomaly pairs in 
each bin. An autocorrelation coefficient matrix will be constructed for the 
combination of different lags. This matrix represents the temporal and spatial 
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variations of the autocorrelation function which is fitted by the Gaussian function 
used in the OTIS model. Finally, the temporal and spatial correlation scales are 
obtained. 

To investigate the vertical variation of the horizontal and temporal scales 
several schemes have been used in deep water analysis. Sprintall and Meyer (1991) 
used the sea surface temperate (SST) and the depth of the 20°C isotherm for the 
Eastern Pacific and White et al. (1982) used 100, 200, and 300 m depths for the 
Northwestern Pacific. Neither scheme is practical for the Yellow Sea region since 
the 20°C isotherm is not always present and the water depth is less than 150 m 
throughout. Since shallow water is strongly influenced by bathymetry, it was 
decided that the seasonal correlations scales for the temperature fields should be 
computed for the surface and depths representing 50% and 80% of the water 
column. 
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n. REGIONAL CHARACTERISTICS 



A. GENERAL 

The Yellow Sea is the middle section of a continental sea system (Figure 2.1). 
To the north is the Bohai Sea and to the south the East China Sea and the Western 
Pacific beyond. It lies between mainland China and the Korean peninsula and covers 
an area of approximately 404,000 km 2 . There are four major rivers that flow into it. 
The Yangtse River located to the southwest, the Yellow and Liao Rivers to the north 
and the Han River to the east. 

The Yellow Sea is considered a shallow basin sea. The bottom slopes from 
20 m in the north to its southern boundary, the shelf break, at 150 m. The Bohai Sea 
region ranges in depth of 20 to 50 m. A trough oriented north to northwest and 
paralleling the Korean coastline occupies the center of the sea. It has a maximum 
depth of 90 m. Within 50 km of the coastline the average water depth is 20 m. 
Moving toward the shelf break the depth slowly increases. The shelf break runs in a 
roughly northeast-southwest orientation and quickly drops from 150 to 500 m in a 
few kilometers. For the remainder of this study the Yellow Sea will include the 
Bohai Sea and that portion of the East China Sea that lies on the continental shelf 
unless specifically stated. 

B. METEROLOGY 

Nester (1977) compiled an extensive report that is the foundation of this 
meteorological overview with wind speeds supplemented by the US Navy Regional 
Climatic Study (Elms, 1989). The Yellow Sea lies within the Asian monsoon 
circulation pattern and experiences hot, short, humid summers and long, cold, dry 
winters. Winter is dominated by an intense cold high pressure system called the 
Siberian high. Its greatest intensity occurs in January, pushing cold, dry air from the 
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Figure 2.1 Bathymetry of the Bohai, Yellow and East China Seas. Depths are 
in 10 m Increments. Bathymetry Data are from the NGDC and 
NAVOCEANO World Wide DBD5 Bathymetry Data Bases. 
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interior of the continent along the east coast of the Asian mainland. The clockwise 
circulation of the Siberian high produces northerly or northwesterly winds referred 
to as the northwest monsoon. This is the predominant feature of the winter weather 
even though small low pressure areas known as Shanghai lows form over the Yellow 
Sea and travel over the Korean peninsula. Figure 2.2 shows the mean pressure and 
wind systems for February. Surface winds are normally from the north-northwest 
with an average speed of 8 m/s. 




Figure 2.2 February Mean Pressure and Surface Wind Pattern During the 
Winter Monsoon (from Langill, 1976). 

March marks the beginning of the transition from winter to spring. The 
weakening of the Siberian high permits the northerly migration of the polar front. 
Tropical maritime air begins to push northward displacing the cold, heavy Siberian 
air causing “migratory highs” that may in turn be displaced by low pressure areas 
from the south. Nort herly winds are no longer dominant and by the end of the season 
winds are from the south averaging 4.6 m/s. 

Summer begins in June and is also referred to as the southwest monsoon 
season. The polar front has migrated to its farthest northern extent. It is a boundary 
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zone between continental air masses from Asia and the warm moist air from the 
tropical oceans and serves as a buffer zone between the two air masses. Figure 2.3 
shows the mean pressure and wind systems for the onset of summer. As summer 




Figure 2.3 June Mean Pressure and Surface Wind Pattern During the 
Summer Monsoon (from Langill, 1976). 

progresses, the Bonin high pressure cell expands westward into the East China Sea 
while the Manchurian low, caused by land mass heating, settles over northern China, 
just north of the Yellow Sea. This is the season of predominantly southward winds 
and heavy rains. Although the three months of summer are ripe for typhoons, the end 
of July and the beginning of September is the period of their greatest probability. 
The East China and central Yellow Seas have a 74% and 32% chance, respectively, 
of a tropical cyclone occurring in the region (Elms, 1 989). The polar front is unstable 
and is conducive to heavy rainfall. The season starts with prevailing winds from the 
south at 5.1 m/s but shift to the north-northeast by September. 

September marks the fall transition season with a rapid southward movement 
of the polar front accompanied by a wide band of rain. By the end of the month the 
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southerly warm, moist, tropical air of the summer monsoon has been replaced by a 
cooler, drier air from the continent. The Siberian high begins to form as the Asian 
continent cools causing cool air masses to periodically break away and be carried 
southward. Winds are predominately from the north at 8.2 m/s. Table 2.1 
summarizes the seasonal cycles of weather phenomena. 

C. OCEANOGRAPHY 

The seasonal changes of the Yellow Sea are strongly influenced by the 
Southeast Asian monsoon. Nester (1977) discusses this monsoonal influence on 
SST, mixed layer depth (MLD) and current speeds and is the basis for this section. 
In winter, when the monsoon winds are at their strongest, the Yellow Sea Current is 
also at its strongest, flowing southward averaging 51 cm/s to the end of the Korean 
Peninsula (Figure 2.4a). Hsueh (1988) found that high northerly wind events during 
winter create a mean sea level gradient that forces bottom water to flow northward 
in spurts along the Yellow Sea trough. Sea surface temperatures in the Bohai Sea 
cool to the freezing point while in the Yellow Sea it is 6°C at its northern extent and 
10 8 C at the shelf break. The MLD is at its deepest during winter owing to both 
convection and wind mixing by the strong monsoon winds. 

March marks the weakening of the northerly winds which in turn leads to the 
weakening of the Yellow Sea Current. The current is variable in direction and speed 
until the end of the season, when it starts flowing northward. The sea surface slope 
is reestablished and a warming trend of the Yellow Sea Current begins (Hsueh, 
1988). The SST’s quickly transition from winter to spring temperatures with an 
average increase of 10°C. The Bohai Sea starts at 6'C and by the end of the season 
it is up to 16°C which is similar to the rest of the Yellow Sea. The MLD is in 
transition so both shallow and deep layers can be found. This is the season of least 
wave height in the Yellow Sea. 
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Polar front 
Soulh of 
R.O.K. 


Heavy 

Rain" 




Rlncking due to slow-movlog 
Mongolian high results In cul-off 
low over Yellnw Sea. 


Pressure 
higher to 
east. 


High pressure area to south; 
low to north. 



Transition 

period 



Under influence of migratory high, 
elongated east -west mean track along 



Heavy snow may he 
exp-cted along Tachaek 
MU with strong north 
to south pressure 
gridient (high to north) 



High pressure ares to west: 
low to east 



Korea under influence of migratory high 4 



Grist winds' 



Tslwan Inws 



Hontal thunderstorm* 
duration i-2 hr*, up to 
twice per month.* 



Typhoons may affect R.O.K* 



Alrmass thunderttorm*. 5-10 km. 
diameter; 2-5 per month. 



Exceptionally good visibility over 
long distance within migratory high 
area. 1 * 



Frontal thunder- 
storms. 0-2/month 



"i cold and 4 warm 
day cycle prevails 
nver Korea. 



Shanghai low m*y develop and Taiwan 

afTcctROK. I lows 

frequently affect R O K 



Sea fog. most frequent In July over Inchnn area. Frontal 
and radiation fog may occur occasionally. f 



Radiation and frontal fogs occur frequently, however duration 
I* short. Radiation fog seen Inland. 



Siberian airmass 
dominates Korea 



1 North Chins low patting over North Korea causes foehn (hot. dry south to cssteriy flow) over west and centrsi sectors of It O K 

* Mongolian low passing over southern MsnrhurU csuies yellow winds (strong, dry west wind carrying yellow dusf). 

5 Ijic frost may occur when migratory high predominates nver R O K. 

* Mean track of migratory high is along J J e N. 

* Low nver Shantung Peninsula may cause gusty winds over R.O.K 

* Intense frnntsi thunderatorma in May signal approach of rainy season i 

I Cool summer and drought if South-west Monsoon weak and Okhotak High Intense. 

* Heavy raln.latter part of rainy season. 

* Hot and humid w|ih maritime tropical air Inmision None or as many at two tropical atorma or typhoons from end nf June thrnugh September may afTert n O K ; least likely In 

August owing to Bonin High. * 

Migratory highs move at about 20 Vnota, speed up with northward movement and vice versa. 

II often precedes had weather 



Table 2.1: Climatic Calendar for South Korea (from Langill, 1976). 
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By the time July arrives the Yellow Sea Current has developed into a closed 
cyclonic circulation (Figure 2.4b) which causes an upwelling effect in die central 

(a) (b) 





Figure 2.4 Representation of Surface Currents Off the Korean Peninsula in (a) 
Winter and (b) Summer (from Bartz, 1972). 

region. This circulation occurs above the thermocline, which is about 13 m in 
August, and may reach a speed of 15.5 cm/s. Below this circulation is a cold water 
mass, commonly referred to as the ‘Cold Water Mass of the Yellow Sea’ (CWYS), 
that remains unchanged and motionless throughout the summer (Li and Yuan, 
1992). The SST is nearly uniform throughout with an average temperature of 23°C. 
The MLD is at its shallowest and maybe negligible in the coastal and the Bohai Sea 
regions due to surface heating. The water is strongly stratified. 

October is the beginning of the transition back to winter conditions. The 
southerly winds have weakened letting the sea surface slope reestablish itself again. 
The Yellow Sea warm current becomes weak and variable. The SST steadily 
decreases 20°C from die beginning of the season to the end. The MLD begins to 
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deepen as stronger winds and convective mixing occur, yet the water mass is still 
stratified. 

Liu et al. (1992) presents a loose interpretation of the temporal and spatial 
distribution of water masses through temperature and salinity analyses. The Yellow 
Sea consists of five main water masses: the coastal water, the Yellow Sea surface 
water, the Bohai- Yellow Sea mixed water, the Yellow Sea bottom cold water and 
the Kuroshio water. The coastal water is influenced by the amount of river runoff. 
In the spring and summer, when run-off is the greatest, low salinities (28.4-30.5) and 
high temperatures (27.5 4 C) identify this water with the reverse true in the fall and 
winter. The Yellow Sea surface water encompasses the middle of the sea. It is the 
dominant water mass in the summer due to the intense surface heating. During the 
fall and winter months the water cools and becomes one water mass with the bottom 
water. The Bohai-Yellow Sea mixed water is in the northern region where the two 
seas connect. This water mass is a mixture of the coastal waters and the Yellow Sea 
surface water. Its greatest extent is in fall and winter when it reaches the Korean 
Peninsula. By summer the Yellow Sea surface water has penetrated into the Bohai 
Sea. The Yellow Sea bottom cold water or the CWYS is present in the spring and 
summer. Its mean temperature is 5 to 12 4 C and the mean salinity 31.8-33.8. As 
stated earlier, with the approach of winter the CWYS merges with the surface to 
form one water mass named the Yellow Sea surface water. The Kuroshio water is 
located along the shelf break and is characterized as high salinity and high 
temperature. The mean temperature is 18-29 # C and the mean salinity is 34.2-34.9. 
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in. OPTIMUM THERMAL INTERPOLATION SYSTEM 



The US Navy requires and conducts daily analyses of the ocean thermal 
structure for a wide variety of tasks. These analyses are made for both global and 
regional areas, with resolutions as small as 20 km in some regions. The Navy’s 
Optimum Thermal Interpolation System (OTIS) combines real-time data, 
climatology and predictions from ocean mixed-layer models to represent an 
accurate picture of the ocean thermal structure on global and regional scales 
(Phoebus, 1988; Clancy et al., 1989). OTIS can be executed in either “shallow” or 
“deep” analysis modes. This statement is somewhat deceptive since the term 
shallow in this case refers to a depth of 400 m. The analysis sequence for OTIS is 
the calculation of (1) SST, (2) the mixed layer depth, and the (3) temperature below 
the mixed layer depth. In truly shallow waters such as the Yellow Sea this analysis 
sequence may be inappropriate since isothermal profiles are possible. 

OTIS is designed to receive its climatological data from the Generalized 
Digital Environmental Model (GDEM) which synthesizes a climatology from a 
given data base. Adequate climatologies for the world’s marginal seas are not 
prevalent but are being rapidly developed by the Naval Oceanographic Office. Even 
though OTIS is not a shallow model in the literal sense, it is hoped that with the 
proper parameters OTIS can be modified for use in marginal seas. 

OTIS is based on the Optimum Interpolation (01) data assimilation technique 
of Gandin (1965). 01 combines a climatological temperature at a gridpoint plus a 
weighted combination of observed and model predicted anomalies. Thus the 
analyzed SST at the k th gridpoint is 
N 

< 31 > 

i = 1 



13 



where T k c is the climatological surface temperature at the gridpoint, Tf the predicted 
temperature at the gridpoint produced by a 24-hour forecast, 7)° the observed 
temperature at location i, 7] c the climatological temperature at location i valid at the 
time of observation i, a fa the weight applied to the observational correction, P* the 
weight applied to the prediction correction, and N the number of observations. With 
the minimal mean square error, these weights satisfy the following equation: 

N 

\i = X ( \ + S ij X< ! ) a kj + ( % + S ik XP k> P* (3 ' 2) 
j= 1 

where % is the space-time autocorrelation function between climatological 
anomalies at location i and location j. V is the noise-to-signal ratio for the 
observation, Vis the noise-to-signal ratio for the prediction and 5 (> is the Kronecker 
delta function. The autocorrelation function is represented as 

\ (33) 
where Ax ijy Ay i; and A/ (> are the east-west, north-south and time separations between 

locations. The quantities A k , B k and C k are the inverses of the east-west, north-south 
and time decorrelation scales at different depth levels, k. These parameters are user- 
specified within OTIS and must reflect the oceans’s resolvable variability about 
climatology. The better the correlation scales are calculated, the higher the accuracy 
the thermal field will be. In order to implement OTIS to the Yellow Sea region 

decorrelation scales, A B~ k x and C ~ k x , must be determined. 
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IV. METHODS OF ANALYSIS 



A. MASTER OCEANOGRAPHIC OBSERVATION DATA SET (MOODS) 
OF THE YELLOW SEA 

1. General 

The oceanographic data used for this study were taken from the Master 
Oceanographic Observation Data Set (MOODS), the U.S. Navy’s environmental 
data base. It is a collection of observed vertical profiles of temperature, salinity and 
computed vertical profiles of sound speed. MOODS consists of over five million 
world-wide profiles of oceanographic data and is considered to be the most 
comprehensive collection of temperature and salinity data known to exist. The data 
base has been compiled from military, academic and commercial sources, both 
foreign and domestic. It has a classified and unclassified version. The data were 
derived from various instrumentation types ranging from the Nansen bottles of the 
1920’s to the present day Conductivity-Temperature-Depth (CTD) recorder. The 
MOODS data are not in a gridded form; they are irregularly distributed in space and 
time. 

The original MOODS data base used for this study consisted of 
approximately 50,000 temperature profiles located between 117-128 °E longitude 
and 25-41°N latitude (Figure 4.1a) which encompasses the Yellow Sea and the 
northern East China Sea. The data spans the period from 1929 to 1991. The data 
were screened to restrict observations only to those temperature profiles that were 
located on the continental shelf. This reduced the data set to 35,658 profiles (Figure 
4.1b). Since horizontal scales have seasonal variations, the data base was binned into 
four seasons. The seasons were defined according to the convention of the Naval 
Oceanographic Office for the Yellow Sea (Haeger, personal communication) 
January, February and March constituted winter; April, May and June, spring; July, 
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August and September, summer; October, November and December fall. Table 4.1 
shows the partition of the data; winter has the least profiles and summer the most. 




a> 

■o 

3 



Figure 4.1 (a) Distribution of Entire MOODS Temperature Profile Data Set, (b) 
Profile Distribution After Removal of Data Deeper than Restricted to the Shelf 
Break. 
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Season 


Number of Profiles 


Winter 


7391 


Spring 


9868 


Summer 


9915 


Fall 


8484 



TABLE 4.1: Number of Temperature Profiles in the Data Base for Each Season. 



The number of temperature profiles from year to year fluctuates greatly. The 
period of 1963 to 1977 is found to have the greatest concentration of profiles 
averaging 315 per season. 
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2. Limitations 

The main limitation of the MOODS data is its irregular distribution in time 
and space. Certain areas and events are over sampled while other areas lack enough 
observations to gain any meaningful insight Vertical resolution and data quality are 
also highly variable depending much on instrument type and sampler expertise. 

The sparsity of profiles have a direct impact on the type of analysis conducted 
on the data. A prime example of this is the Bohai Sea region in winter. Access for 
data gathering in this region has been limited. In fact, most observations in the area 
have been conducted prior to the 1960’s. Figure 4.2a depicts unclassified 
observation locations for the winter seasons of 1963 to 1977. Another data poor area 
is the eastern coastline of mainland China (Figure 4.2b). The spring of 1969 had the 
greatest number of profiles for a given season yet the entire Chinese coastline lacks 
a single recorded hydrographic cast. An analysis of location plots from 1963 to 1977 
(Appendix A) shows this same general lack of data for these two regions. Statistical 
relevance must be considered when trying to calculate spatial scales with such 
sparsely located observations. 

Temporally uneven distributions are just as great. Over all, the data base 
ranges from a low of one profile for 1944 to a high of 1790 profiles in 1968. The 
number of yearly temperature profiles in the 15 year analysis period (1929 to 1991) 
is ranges from 592 to 925. When further subdivided into seasons, this greatly 
reduces the pool of temperature profiles for which statistical analysis may be 
conducted. Figure 4.3 and Appendix B provide examples of the temporal 
irregularities and sparseness of the observations. There are frequent 15 to 20 day 
gaps of ten observations or less per day. For an area as large as the Yellow Sea this 
does not provide an adequate data base to establish the inherent variability on a 
yearly basis. 



17 



(a) 



(b) 





Longitude (E) 

Figure 4.2 (a) Winter Temperature Profile Locations from MOODS Data Base From 
1963 to 1977, (b) Locations of Temperature Profiles During Spring 1969, a Period 
Illustrating Maximum Data Density in a Given Year. 




Figure 4.3 Number of Temperature Observations During Winter and Spiing 1977, 
Illustrative of Typical Data Density in Any Given Half Year. 
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Spatial and temporal irregularities along with the shear lack of data must be 
carefully weighed in order to avoid statistically-induced variability. Originally, it 
was hoped that there would be sufficient data to compute horizontal autocorrelation 
scales on a yearly basis and determine a trend. After analyzing the location and 
histogram plots (Appendix A and B), it is obvious that statistically significant results 
can not be obtained on a yearly basis. Analysis was done on a mean seasonal basis 
using the data for all years. 

B. CLIMATOLOGY OF THE YELLOW SEA TEMPERATURE FIELD 

1. General 

Building a horizontal temperature field climatology from the MOODS data 
base is not a straightforward task and in shallow water the problems are magnified. 
Most of the problems relate to lack of data in contrast to temporal and spatial 
variability. Without enough data single observations force artificial gradients to 
appear. Seasonal time scales are thought to be more important than spatial scales but 
shallow water mitigates this by introducing much smaller temporal scales. In 
addition, the sloping bottom of the Yellow Sea poses difficulties with interpolating 
techniques. 

2. Method for Establishing Temperature Climatology 

Four main steps were used in establishing the temperature climatology: (1) 
seasonally binning the data, (2) interpolating the temperature profiles for the desired 
depth, (3) gridding the temperature with the Minimum Curvature Spline Interpolator 
(MISP) algorithm and (4) determining the desired grid spacing. The entire MOODS 
data base was binned by three month seasons as described earlier and then the 
temperature profiles were interpolated for the specified depths. Climatology 
building in deep water is distinctly different to that of shallow water. To compute 
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deep water climatologies several submodels are constructed in the vertical (Figure 
4.4a). These submodels extend temperature profiles that do not extend to the sea 
floor (or to a given depth level) by interpolating values from profiles that do extend 
each submodel level. 

In shallow waters having sloping bottoms, if vertical submodels are used 
(Figure 4.4b), some areas (hashed area) are not covered by any submodels and 
therefore no climatology can be constructed for the area (Haeger, 1991). In order to 
compensate for this problem a vertical interpolating scheme that produces 
horizontal sloping surfaces was used (Figure 4.5). Horizontal temperature 
climatologies were constructed for the surface, 50% and 80% depth of the water 
column. Mid-water properties were expected to be modeled by the 50% depth 
whereas, the 80% depth was chosen since it was near the bottom but provided 
enough height for small scale depth variations. The surface temperatures in this case 
were actually any observations from the surface to 2 m below. 

After the temperature profiles were interpolated for the specified depths the 
temperature values were then horizontally gridded using the Minimum Curvature 
Spline Interpolator (MISP) algorithm. The MISP can be broken into three stages: 
partial grid computation, regional grid computation and final computation. The 
partial grid computation consists of binning the input data into bins equal to the 
specified output grid size. All data falling within a distance of half a grid cell of a 

grid point (i, j) were averaged to form an averaged x.j, ^.(location) and 7V. 
(temperature) set associated with that grid cell. At this stage 7V is not situated at the 
gridpoint (i, j ) but is located at the position (x^ y.J . After all data have been read, 
a ‘partial grid’ has been formed and consists of grid points with and without 
associated ( 7V, x.j, y.p values. 
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The second stage of the gridding process consists of computing a regional 



grid. The (7V, * , y. tJ ) sets obtained for a grid cell in the first stage are used as input. 




Figure 4.4 The Geometries of Vertical Submodel in the Construction of Climatology 
Temperature Profiles in the (a) Deep Water and (b) Continental Shelf. 




surface 



50% 

80% 



Figure 4.5 Depiction of the Depth Levels that Horizontal Temperature 
Climatologies were Constructed for in the Yellow Sea. 

These values are first smoothed by rebinning to a grid three times the size of the 
original grid. A grid cell may or may not have a value at this stage. The minimum 

curvature spline interpolation routine was used to interpolate temperature (TV) 
values to the gridpoints (i, j ) centered at the grid cell resulting in the completion of 
the regional grid. 
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The third stage is to adjust the mean temperature ( 7V) values of the partial 

grid to their respective gridpoints (/, j). First, the regional grid is interpolated back 
to the original grid size with a bi-cubic spline to match the partial grid. Then the 
gradient computed for each regional grid cell to the associated partial grid cell. The 
new partial grid is then merged with the regional grid. The merging consists of using 
the partial adjusted grid values where they exist at a gridpoint and the regional grid 
values where the partial grid cells are empty. This allows the input data to dominate 
in data rich areas and the regional data to dominate in data sparse areas (Haeger, 
personal communication). 

Sometimes a climatology field is sensitive to the selection of the grid size. 
Too small a grid spacing produces false gradients appearing as bull’s-eyes. Yet, if 
the spacing is too large, features that do exist may not appear. Grid spacings of 20 
minutes and 30 minutes were analyzed in this study. As expected, the contours 
associated with a 30 minute grid are smoother (Figure 4.6a) than the 20 minute grid 
(Figure 4.6b) causing mesoscale features to be lost. In general, the 20 minute grid 
better depicts the findings of Elms (1989) for the Yellow Sea and hence this grid will 
be used for the rest of the study. 

3. Yellow Sea Mean Temperature Fields 

The surface temperatures for all four seasons (Appendix C) show that 
temperatures increase toward the south. During winter, freezing temperatures occur 
in the Bohai Sea region which result in surface ice. In the Yellow Sea winter 
temperatures are significantly higher, with readings in the high teens (14-18 °C) near 
the shelf break. As the year progresses towards summer, the horizontal temperature 
gradient decreases. Spring temperatures range from 6 °C in the Bohai Sea, 14 °C in 
the central Yellow Sea and 24 °C near the Kuroshio Current. By the time summer 
arrives the SST of the entire Yellow Sea is above 20 °C and is almost uniformly the 
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same. With the onset of cooler fall conditions the north-south gradient reappears. 
The Bohai Sea shows the presence of a pool of warm (18 °C) water while the Yellow 
Sea begins to establish its winter gradient with lower temperatures (14 °C) to the 
north and a steady increase to the shelf break (24 °C). 



(a) 




110 120 122 124 126 128 



Longitude ( E) 



(b) 




Figure 4.6 Climatological SST Field on (a) a 30 Minute Grid and (b) a 20 Minute 
Grid. 



The horizontal temperature contours at 50% and 80% water column depths 
tend to mimic each other from season to season. The winter subsurface 
climatologies closely mirror that of the surface. The largest difference is in the shelf 
break area where surface temperatures are slightly higher. On the other extreme, 
summer shows the most contrast between surface and subsurface temperatures. The 
surface temperature field shows almost uniform values throughout whereas the 50% 
and 80% depth contours show a cold core in the middle of the Yellow Sea. Both 50% 
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and 80% depth levels show this cold core beginning to form in spring and 
dissipating in fall. 

C. SPACE /TIME SORTING 

At present, the spatial and temporal scales for the OTIS model are artificially 
prescribed as 100-250 km (Phoebus, 1988) and 20 days (Clancy et al., 1989) for the 
region of Japan and the South China Sea. Since synoptic temporal and horizontal 
scales are usually smaller than 30 days and 300 km (Mooers and Robinson, 1981), 
it was decided that a matrix with a maximum temporal lag of 30 days with a one day 
increment and the maximum spatial lag of 300 km with a 10 km increment would 
be used. Each matrix bin is associated with a specific temporal and spatial lag. Other 
matrix sizes were tested. It was found that smaller temporal and spatial lags 
excluded too much data while larger lags did not necessarily improve the results. 

As mentioned earlier, an important feature of the Yellow Sea thermal field is 
the existence of the CWYS (Li and Yuan, 1992). This water mass is located in the 
central bottom of the Yellow Sea. In the winter it merges with the surface water but 
from April to November it is a distinct water mass with temperatures of 4-5 °C. 
During the summer the overlying waters are 10-15°C warmer. The isotherms are 
closed curves (Figure 4.7), suggesting the basin may be treated as an isolated system 
(Guan, 1994). Therefore, it is reasonable to assume that the spatial autocorrelations 
of temperatures inside the Yellow Sea basin are isotropic, i.e.,the autocorrelations 
are independent of direction (N-S or E-W) and depend only on the distance between 
two temperature locations. This assumption of horizontal isotropicity greatly 
simplifies the computation. 

After the matrix size was defined and the isotropic assumption made, 
anomaly pairs were placed within the temporal and spatial bins. To accomplish this, 
it is necessary to take each observation,^, find the closest gridpoint and compute 
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the anomaly, T o \ by subtracting the climatological value at this gridpoint T k from 
the observation, i.e., 7* ' = T - T k . For every individual temperature anomaly, T o \ 
it is paired with the other data points, T 0 \ within the data base. The temporal and 




Figure 4.7 Temperature Distribution in Summer (April-November) of the ‘Cold 
Water Mass of the Yellow Sea’ (CWYS) Near the Sea Floor (from Su and Weng, 
1994). Isotherms in Degrees Celsius. 



spatial difference or lag is then calculated between the two anomalies. The anomaly 
pairs (T o \ t a ') are then placed in the corresponding temporal and spatial lag bin. 
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Three separate data runs were conducted. The first data run was to use 
separate seasons and individual years from 1963 to 1977. The second run was to use 
all the MOODS data binned by season. There was no delineation between years; 
therefore an observation from 23 January 1963 could be paired with an observation 
from 24 January 1969 and only have a lag time of one day. The third run also 
included the entire data set binned seasonally but each observation retained its 
year.Hence, it was paired with observations of the same year, then sorted into the 
matrix. For simplicity the different input types will be referred to as Type I, Type II, 
and Type HI. 

The number of observations in the data set (AO will determine the number of 
pairs that can be binned. The number of pairs ( S ) is computed by 

5 = ay-n (4.D 

2 

In Type I and Type II data the computation S is straight forward but for Type HI the 
value of S is the summation of S for each individual year. As to be expected, Type I 
has the fewest pairs in each bin. It ranges from one pair to a high of 2000 pairs. In 
contrast Type II ranges from approximately 750 to over 24000 pairs and Type HI 
from approximately 40 to the 4000. Figure 4.8 shows the pair distribution of Type I 
and Type H. It can be easily seen that the Type I scheme (Figure 4.8a) produced 
three main (green) concentrations of profile pairs, the other areas lack substantial 
data. Type H (Figure 4.8b) has a much higher and more even distribution. From this 
lack of data in Type I it was determined that statistically relevant information could 
not be obtained. This study will analyze the advantages and disadvantages of Type 
II and Type III data. 
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Figure 4.8 Contour Plots of Time-Distance Lag Pair Distribution (a) TYPE I for 
Spring of 1969 and (b) TYPE II for Spring. 
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D. AUTOCORRELATION FUNCTION 



After the anomaly pairs have been spatially and temporal sorted, an 
autocorrelation function (r| (m,n)) for each bin with temporal lag nAt and 
horizontal lag mAy (At = 1 day, Ay = 10 km) identified as the bin (m, n) is computed 
by 



H (m, n) = 



y T 'Tn 

jLu o o 

bin (m, n) 



/2 



,2 



(4.2) 



o 



Hbin(m,n ) *ibin(m,n) 



Figure 4.9 shows the temporal variation of autocorrelation ( ( m , n) r, at the 
surface for the winter season of Type II and Type HI data. The three curves on each 
graph show the dependency of r\ (m, n) on temporal lag nAt for three different 
spatial lags: m = 0 (no lag, dash-dot), m = 1 (10 km lag, solid), and m = 15 (150 km 



lag, dots). Such a decrease of r\ (m, n) with nAt reveals approximately a Gaussian 
type relation. On both graphs there is no unity or a correlation of one at zero lag. This 
occurs since the covariance estimates at zero temporal lag did not include products 
of each data point with itself and also because substantial subgrid noise was present 

(White et al., 1982). r\ (0, 0) is greater for Type HI than the Type II computations 
which is expected since temperatures of the same year are being compared with each 
other. Type II plots (Figure 4.9a) are much smoother than the Type HI plots (Figure 
4.9b). This high frequency temporal variability may be induced by a smaller data 
base and by forcing mechanisms like meteorological, tidal and freshwater runoff 
that were averaged out in the Type II computations but are clearly recorded in the 
Type III calculations. It is unclear how many or what type of natural forcing 
mechanisms are creating this variability or if it is an artifact of having too little data 
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that artificially induces variability. In each case the e-folding time for zero lag 
distance is approximately 13 days. 

The spatial autocorrelations showed similar characteristics (Figure 4.10). 
Type III had higher autocorrelations at zero lag. Variability is large in the Type HI 
scheme and the reason for this is still unclear. The e-folding distance for Type II data 
was roughly estimated from Figure 4.10 to be approximately 125 km while for the 
Type III it is 150 km. Since the e-folding scales are similar for both computations 
(Type II and Type IH) and Type II exhibits statistical stability, the Type II 
computations were considered the more appropriate choice in calculating the 
decorrelation scales. 

E. E-FOLDING DECORRELATION SCALES OF THE YELLOW SEA 
THERMAL FIELD 

With regard to the surface and subsurface thermal indicator parameters, the 
e-folding temporal and spatial decorrelation scales can be obtained from the 
autocorrelation coefficient function, r| (m, n ) . The variations of ri (m, n) verses the 
horizontal spatial and temporal lags (mAy and nAt) reveals a Gaussian-type 

dependency (e.g., Figure 4.9). The autocorrelation matrices T| (m, n ) are fitted by 

the Gaussian function 

-A 2 (mAY) 2 -Cj(„A<) 2 

T 1 W (m,n) = 11(0,0)* (4,3) 

where k denotes the different levels (k =1 for the surface, k =2 for the 50% water 
column depth, and k =3 for the 80% water column depth). The regression 

coefficients, A~ l and C^ 1 , are the e-folding horizontal spatial and temporal 

decorrelation scales. For a given time lag nAt, the coefficient A^ 1 is interpreted as 
an e-folding distance, m i Ay, for the level k 
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(a) Winter 




Time lag (days) 



(b) Winter 




Figure 4.9 Temporal Variations of Autocorrelation Function r\ (m, n) for (a) Type 
II and (b) Type III Data at the Surface.The Three Curves on Each Graph are the 
Dependency of t\ (m, n) on Time Lag nAf for m = 0 (no lag, dash-dot), m = 1 (10 km 
lag, solid), and m- 15 (150 km lag, dots). Data were Obtained from the MOODS 
Temperature Data Base for the Years 1929 to 1991. 



30 



(a) Winter 




(b) Winter 




Figure 4.10 Spatial Variations of Autocorrelation Function t| (m, n ) for (a) Type II 
and (b) Type III at the Surface. The Three Curves on Each Graph are the 
Dependency of T| (m, n ) on Spatial Lag (m Ay) for n = 0 (no temporal lag, dash-dot), 
n = 1 (1 day lag, solid), and n = 15 (15 day lag, dots). Data were Obtained from the 
MOODS Temperature Data Base for the Years 1929 to 1991. 
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(4.4) 



T| (m, n) 

-A 2 (m.AY) = In— = -1 

T| W (0, n) 



(k) 

where the autocorrelation function reduces from zero spatial lag, r\ (0, n ) , to 
m.A Y spatial lag, T| W (m,., n) = ^ _ 1 T| w (0, n) . 

Similarly, for a given spatial lag (mAy ) , the coefficient C~ k is interpreted as 
an e-folding temporal separation (nAt) for the level k 



2 2 “H ( m > n ;> 

-w {nAt) 2 = In — — - 

* T| ^ (m, 0) 



= -1 



(4.5) 



( k) 

where the autocorrelation function reduces from rj (m, 0) to 
rj (A:) (m, n f ) = e~ l r \ W (m, 0) . 



F. SIGNAL-TO-NOISE RATIO 

In applying the optimum interpolation for the OTIS model or an optimal 
network design, certain statistics of the field of interest such as the variance of both 
the signal and the noise must be known. The signal variance for the bin (m, n) is 
estimated by 

a 2 , x = Y T'l (4.4) 

s(m,n) Au o v 

bin ( m,n ) 

The noise variance is derived from two sources (Sprintall and Meyers, 1991), 
geophysical error and instrument error. Geophysical error is unresolved thermal 
variability with scales smaller than the typical time and distance scales between two 
temperature profiles. Instrumental error can be shown to be negligible compared to 
geophysical noise (White and Bernstein, 1979). These noise sources inhibit the 
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detection of large-scale variability and can be an important factor in establishing the 
sampling density necessary to detect the signal (Sprintall and Meyers, 1991). 

In calculating the autocorrelation value of the first bin (zero to 1 day and zero 
to 10 km lag) the temperature anomalies were never paired with themselves. The 
values of the autocorrelation coefficient function at zero temporal and spatial lag 

[t| (0, 0) ] were extrapolated from the first bin into the origin. Therefore, r| (0, 0) is 
not equal to 1.0. The relationship between the temperature signal and the noise 
signal can be described by the signal-to-noise parameter (Sprintall and 
Meyers, 1991) 



X = 



3 (Q, Q) 

1-ti(0,0) 




(4.5) 



where o 2 is the noise variance and a 2 is the signal variance. X represents the signal- 
to-noise ratio or level. The larger the X, the less noise that exists in the bin. 
Theoretically, if there is no noise, r| (0, 0) = 1 and o 2 = 0. Therefore X will be °°. 



O 0 

If X > 2, the ratio of the signal variance, o s , to the noise variance, o n , is greater than 
4 which is considered good by White et al. (1982) and Sprintall and Meyers (1991). 
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V. DISCUSSION 



A. TEMPORAL AND SPATIAL DECORRELATION SCALES 

All temperature profiles from 1929 to 1991 were used in the calculation of 
temporal and spatial decorrelation scales. Delineation of separate years within the 
autocorrelation computations was not permitted (Type II averaging) which averaged 
the various forcing mechanisms prevalent in this shallow water basin. The computed 
autocorrelations for each season were then fitted to a Gaussian function of the form 
of (4.3). The regression method was used to calculate the spatial and temporal 

decorrelation scales at three different levels, A k and C^ 1 , respectively. Since an 

isotropic assumption concerning the spatial temperature distribution was made, 
there was no directional dependency. This means that decorrelation scales were 
assumed to be the same throughout the Yellow Sea basin. Table 5.1 presents these 
decorrelation scales and signal-to-noise ratios for winter and summer. 



Season 


Depth 


Temporal 

(days) 


Spatial 

(km) 


Signal-to 
noise ratio(7i) 


Winter 


sfc 


14.7 


158 


2.06 




50% 


14.7 


167 


1.88 




80% 


15.2 


172 


2.21 


Summer 


sfc 


12.3 


251 


3.0 




50% 


15.8 


169 


2.21 




80% 


17.2 


157 


3.0 



Table 5.1: Seasonal Temporal and Spatial Scales 



The signal-to-noise ratio for both winter and summer are higher than the 
signal-to-noise ratio for the deep waters of the Eastern Equatorial Pacific (~ 1.0) as 
reported by Sprintall and Meyers (1991). This infers that the temperature signal in 
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this study is stronger (more stable) than the background noise. In both seasons the 
smallest signal occurs at the 50% water column depth. This is expected since this is 
the transitional (thermocline) depth between the surface layer and the deep layer, 
where meteorological and topographic effects occur. The surface signal-to-noise 
ratio was smaller in winter than in summer. 

Winter temporal and horizontal decorrelation scales are the smaller than 
summer and exhibit nearly uniform values at all three depths. The temporal scales 
vary only by a half a day and the horizontal scales by only 14 km. This uniformity 
is expected considering the oceanographic conditions of the region. Sea surface 
icing occurs to the north in the Bohai Sea and the northern Yellow Sea and along the 
northern coastal regions of Korea and China. Temperatures range from freezing in 
the north to 20°C in the south along the shelf break. This large temperature range 
will inherently make horizontal decorrelation scales small. Vertically, winter 
experiences the greatest mixing (Nester,1977) with monsoonal winds at their 
strongest throughout the entire season. With a large part of the Yellow Sea having 
depths less than 50 m, mixing reaches to the bottom creating isothermal profiles. 
Thus, decorrelation scales will be similar from the surface to the 80% water column 
depth. 

Summer surface spatial decorrelation scales are 100 km longer than winter. 
This is reasonable since this region experiences large solar heating in the summer. 
The summer climatology (Appendix C, Figure C.l) shows the SST to be nearly 
constant throughout the region. The surface temporal scale is less by five days than 
at deeper depths. This may be due to a combination of factors. Besides the large 
solar heating that occurs in the summer, the summer monsoonal winds are much 
weaker than the winter winds. Therefore, vertical mixing is not as great as in winter. 
A strong thermocline is present in summer. Thus, the upper layer water is affected 
by the meteorological synoptic events rather than the entire water column as is in 
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winter. The temporal scale at the surface should be shorter in summer than in winter 
(Table 5.1). In the deeper layers (50% and 80% depths) temperatures are not directly 
affected by the meteorological events of the surface and experience longer temporal 
scales. The 80% water column depth has the longest temporal scale which could be 
directly related to the ‘Cold Water Mass of the Yellow Sea’ (CWYS) (Li and Yuan, 
1992) that exists throughout the summer in the middle of the Yellow Sea. Spatial 
scales are close to 100 km smaller than that at the surface. Without the extreme 
vertically mixing, the CWYS retains some its winter water mass characteristics from 
which it is derived (Liu et al., 1992) and remains a distinct cold water mass that 
resists mixing. This will cause smaller decorrelation scales. 

B. OBSERVATIONAL NETWORK DESIGN 

The temporal and horizontal decorrelation scales provide useful information 
for an optimal design of an observational network. The principle is to minimize both 
instrumental and geophysical errors. Since noise due to instrumental error in XBT 
measurements is usually only about O.l’C (Barnett and Bernstein, 1980), and even 
smaller in CTD measurements, instrumental error can be neglected. This implies 
that the reduction of noise must be accomplished by the reduction of geophysical 
error. This is usually done by increasing the sample density. Sprintall and Meyers 
(1991) suggested that an optimal sampling density of two or three samples per 
decorrelation scale would be adequate to resolve thermal variability. 

This would mean that spatially, any future temperature measurement in the 
winter may be conducted at 50-80 km intervals with the knowledge that the 
subsurface features will also be adequately sampled. The summer season will 
require more consideration. If only surface features are considered, a larger 
measurement interval of 80-125 km may be used. In order to provide full coverage 
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of subsurface features a 50-80 km interval is recommended.Temporally, a 4-6 day 
sampling interval will provide the min imal coverage for both seasons. 

C. UNCERTAINTIES IN THE TRANSITIONAL SEASONS 

The transitional seasons, spring and fall, showed unusual temporal and 
spatial decorrelation scales (Table 5.2). The scales were significantly larger than in 
winter or summer. The spatial surface scales computed are 521 km in the spring and 
757 km in the fall. It might be expected that the transitional surface scales be longer 
than winter and summer due to the lack of meteorological synoptic forcing or the 
tendency for the temperature contours (Appendix C), especially in spring and to a 
lesser degree in fall, to rap around or follow the coastline around the Yellow Sea 
basin. This could mean that an eastern coastal temperature profile would be more 
correlated to a western coastal profile 300 km away than it is to the profiles in the 



Season 


Depth 


Temporal 

(days) 


Spatial 

(km) 


Signal-to 
noise ratio(X) 


Spring 


sfc 


20.5 


521 


4.9 




50% 


17.7 


226 


2.48 




80% 


18.0 


204 


2.38 


FaU 


sfc 


18.0 


757 


4.9 




50% 


16.9 


278 


2.48 




80% 


18.2 


179 


2.59 



Table 5.2: Seasonal Temporal and Spatial Scales 



middle of the basin that are closer. In either case, the extremely large scales are 
unrealistic if not for the only reason that the enclosed Yellow Sea basin itself is 
approximately 450 km in the east-west direction and over 700 km in the north-south 
direction. 
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The temporal scales are approximately six days longer than in winter and 
summer. There was no evidence to support these larger values. Since the onset or 
cessation of seasonal changes vary from year to year and this study specifically 
defined April, May and June as spring and October, November, and December as 
fall, statistical errors are expected. It is highly recommended that future studies 
examine different month groupings for these transition seasons. 

Vertically both the temporal and spatial scales decreased with depth. The 
larger change is in the spatial scales. In spring there was a 300 km and in fall a 500 
km difference from the 50% and 80% depths to that at the surface. As stated earlier, 
these findings seem highly improbably and the irregularities of the MOODS data 
and the determination of season must be examined further. 

D. LIMITATIONS 

The purpose of this study was to calculate the seasonal temporal and spatial 
decorrelation scales to be used by OTIS. Scales were obtained but the validity of 
these findings must be examined and a potential user discouraged from blindly 
implementing these scales into various models. 

The isotropic assumption of uniform conditions throughout the basin was 
based on the existence of closed isotherms for the majority of the time in the Yellow 
Sea basin. This made computations feasible and provided a general reference for 
temporal and horizontal scales. In reality, there probably exists a similar latitudinal 
and longitudinal dependence as in deep water (White et al., 1982; Sprintall and 
Meyers, 1991). This is due to the fact that isotherms are not closed year round, the 
basin is elongated in the north-south direction and coastal regions typically have 
different cross-shore and along-shore scales. Thus, future studies should examine 
the possibilities of calculating two directional scales. 
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The convention used in defining seasons can have a large statistical impact. 
The convention used in this study was recommended by the Naval Oceanography 
Office. It was found that spring and fall had unrealistic surface scales which may be 
a direct result of the seasonal binning. A more realistic approach may be to define 
spring and fall as two months or less in duration as opposed to three used in this 
study. 

The MOODS data base is a limiting factor of itself. In this study, 
decorrelation scales that pertain to the entire Yellow Sea basin were desired. Yet, the 
data were biased in location (Appendix B). Most of the profiles were located in the 
eastern region of the sea, with significantly fewer in the Bohai Sea and along the 
coastline of China. Temporally the profiles were irregular (Appendix A) from 
month to month and year to year. This irregularity led to the analysis of three 
different computational schemes. 

The first scheme, Type I, calculated autocorrelations on a yearly basis. 
Values were erroneous since not enough data existed to produce statistically 
relevant results. Type II calculations used the entire MOODS data base but no 
delineation was permitted between years. This scheme averaged the high frequency 
variability such as meteorological events, tidal forcing and river runoff. This scheme 
showed a definite seasonality (Appendix D) both temporally and spatially. Type ID 
calculations also used the entire MOODS data but this time temperatures were 
paired with only temperatures of the same year. Although the autocorrelations 
computed from Type II and Type III schemes are similar, the Type IE scheme 
produced large variability both temporally and spatially (Appendix D). It is not 
known if the variability is statistically produced by a smaller data base or are the 
result of the many shallow water forcing mechanisms. Therefore, more effort is 
needed to determine the best use for both schemes. 
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The last limitation to be discussed is the fitting function. Even though it is 
customary to fit the calculated autocorrelations to a Gaussian curve, as in the OTIS 
model, it may not be the most accurate procedure for shallow water. The winter and 
summer autocorrelations rapidly decreased in time and space. The spring and 
summer curves did not. In future studies it may be better to use a combination of a 
linear function and a Gaussian curve. 
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VI. CONCLUSION 



The goal of this study was to compute seasonal temporal and spatial 
temperature autocorrelation scales for the Yellow Sea. These scales are an integral 
part of the accurate running of the OTIS model which at this time is designed for 
deep water thermal analysis. Decorrelation scales were calculated for three depths, 
i.e., the surface and at 50% and 80% of the water column depth, and for four seasons. 
Variations were found to be both seasonally and depth dependent. 

The temperature data were comprised of approximately 36,000 observed 
profiles from the MOODS worldwide data base. All profiles were located on the 
continental shelf and spanned the time period of 1929 to 1991. The widely irregular 
distribution of temperature profiles in time and space made yearly analysis 
statistically irrelevant. Two time-organized schemes were analyzed: (1) 
temperatures were paired with no delineation between years, referred to as Type II 
data and, (2) temperatures paired with temperatures of the same year, referred to as 
Type HI data. The Type II scheme was used in the calculation of the temporal and 
spatial scales since it produced a larger data set and averaged the high frequency 
variability due to individual forcing mechanisms. 

The winter season had the smallest temporal and spatial scales of the four 
seasons. This is attributed to the strong cold monsoonal winds which produced the 
strongest vertical mixing of all the seasons and the extreme temperature range 
extending from freezing in the north to 20°C in the south along the continental shelf. 
The temporal decorrelation scales approximated to 15 days and spatial decorrelation 
scales to 165 km at all depths. 

The summer season had the next smallest scales, but unlike winter it showed 
some vertical variability. The scales range from 12.3 days and 251 km at the surface 
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to 17.2 days and 157 km at 80% depth. The approximately 100 km difference 
between the levels is due to the lack of vertical mixing and the presence of CWYS. 

An optimal network design was suggested. It was found that a 50-80 km 
interval will provide adequate sample density to discern winter features. The 
summer season exhibits a surface and subsurface minimum density interval. To 
discern only surface features a 80-125 km interval is adequate but for complete 
depth coverage a 50-80 km interval should be used. Temporally, a 4-6 day sampling 
interval will provide coverage for both seasons. 

The use of these scales in the OTIS model must be cautiously considered. The 
inherent problems in the MOODS data base, irregularity in space and time, the 
organization of the data base into seasons and the neglect of different forcing 
mechanisms all have an effect on the calculation of the temporal and spatial scales. 
More temperature profiles and analysis is required in order to better refine these 
scales and improve the OTIS model. 
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APPENDIX A SEASONAL LOCATION PLOTS OF MOODS 
TEMPERATURE PROFILES FROM 1963 TO 1977 
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Figure A.l Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
for 1963. 
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Figure A.2 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
for 1964. 
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Figure A.3 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
for 1965. 
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Figure A.4 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
for 1966. 
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Figure A.5 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
for 1967. 
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Figure A.6 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
for 1968. 
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Figure A.7 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
for 1969. 
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Figure A.8 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
for 1970. 
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Figure A.9 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
for 1971. 
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Figure A. 10 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
for 1972. 
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Figure A.l 1 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
for 1973. 
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Figure A. 12 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
for 1974. 
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Figure A. 13 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
for 1975. 
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Figure A.14 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
for 1976. 
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Figure A. 15 Seasonal Location of MOODS Temperature Profiles in the Yellow Sea 
for 1977. 
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APPENDIX B TEMPORAL DISTRIBUTION OF MOODS 
TEMPERATURE PROFILES FROM 1963 TO 1977 
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Figure B.l The Number of MOODS Temperature Profiles During 1963 in the 
Yellow Sea. 
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Figure B.2 The Number of MOODS Temperature Profiles During 1964 in the 
Yellow Sea. 
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Figure B.3 The Number of MOODS Temperature Profiles During 1965 in the 
Yellow Sea. 
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Figure B.4 The Number of MOODS Temperature Profiles During 1966 in the 
Yellow Sea. 
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Figure B.5 The Number of MOODS Temperature Profiles During 1967 in the 
Yellow Sea. 
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Figure B.6 The Number of MOODS Temperature Profiles During 1968 in the 
Yellow Sea. 
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Figure B.7 The Number of MOODS Temperature Profiles During 1969 in the 
Yellow Sea. 
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Figure B.8 The Number of MOODS Temperature Profiles During 1970 in the 
Yellow Sea. 
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Figure B.9 The Number of MOODS Temperature Profiles During 1971 in the 
Yellow Sea. 
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Figure B.10 The Number of MOODS Temperature Profiles During 1972 in the 
Yellow Sea. 
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Figure B.ll The Number of MOODS Temperature Profiles During 1973 in the 
Yellow Sea. 
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Figure B.12 The Number of MOODS Temperature Profiles During 1974 in the 
Yellow Sea. 
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Figure B.13 The Number of MOODS Temperature Profiles During 1975 in the 
Yellow Sea. 
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Figure B.14 The Number of MOODS Temperature Profiles During 1976 in the 
Yellow Sea. 
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Figure B.15 The Number of MOODS Temperature Profiles During 1977 in the 
Yellow Sea. 
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APPENDIX C YELLOW SEA SEASONAL MEAN TEMPERATURE 

FIELDS 
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Figure C.l Yellow Sea Seasonal Mean Temperature Field for the Surface on a 20 
Minute Grid. Data were Obtained from the MOODS Temperature Data Base for the 
Years 1929 to 1991. 
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Figure C.2 Yellow Sea Seasonal Mean Temperature Field for 50% Water Depth on 
a 20 Minute Grid. Data were Obtained from the MOODS Temperature Data Base 
for the Years 1929 to 1991. 
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Figure C.3 Yellow Sea Seasonal Mean Temperature Field for 80% Water Depth on 
a 20 Minute Grid. Data were Obtained from the MOODS Temperature Data Base 
for the Years 1929 to 1991. 
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APPENDIX D AUTOCORRELATION FUNCTIONS 
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Figure D.l Temporal Variations of Autocorrelation Function ri (m, ri) for Type II 
Data at the Surface.The Three Curves on Each Graph are the Dependency of 
T| (m, ri) on Time Lag nAf form = 0 (no lag, dash-dot), m = 1 (10 km lag, solid), and 
m = 15 (150 km lag, dots). Data were Obtained from the MOODS Temperature Data 
Base for the Years 1929 to 1991. 



85 



Winter 



Spring 



0.8 




-0.2 



0 10 20 30 



Summer 





o 



- 0.2 



0 10 20 30 



Fall 




Time lag (days) 



Figure D.2 Temporal Variations of Autocorrelation Function r\ ( m , n) for Type II 
Data at 50% Water Depth. The Three Curves on Each Graph are the Dependency of 
ri (m, n) on Time Lag nAr for m = 0 (no lag, dash-dot), m = 1 (10 km lag, solid), and 
m = 15 (150 km lag, dots). Data were Obtained from the MOODS Temperature Data 
Base for the Years 1929 to 1991. 
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Figure D.3 Temporal Variations of Autocorrelation Function T| (m, n) for Type II 
Data at 80% Water Depth. The Three Curves on Each Graph are the Dependency of 
t | (m, n) on Time Lag nAr for m = 0 (no lag, dash-dot), m = 1 (10 km lag, solid), and 
m = 15 (150 km lag, dots). Data were Obtained from the MOODS Temperature Data 
Base for the Years 1929 to 1991. 
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Figure D.4 Temporal Variations of Autocorrelation Function Tj (m, n ) for Type HI 
Data at the Surface. The Three Curves on Each Graph are the Dependency of 
T| (m, n ) on Time Lag nA/ for m = 0 (no lag, dash-dot), m = 1 (10 km lag, solid), and 
m = 15 (150 km lag, dots). Data were Obtained from the MOODS Temperature Data 
Base for the Years 1929 to 1991. 



88 



Winter 



Spring 




Summer 




Time lag (days) 




Fall 




Time lag (days) 



Figure D.5 Temporal Variations of Autocorrelation Function r| (m, n) for Type HI 
Data at 50% Water Depth. The Three Curves on Each Graph are the Dependency of 
H (m, n) on Time Lag nAf for m = 0 (no lag, dash-dot), m = 1 (10 km lag, solid), and 
m = 15 (150 km lag, dots). Data were Obtained from the MOODS Temperature Data 
Base for the Years 1929 to 1991. 
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Figure D.6 Temporal Variations of Autocorrelation Function T| (m, n) for Type IQ 
Data at 80% Water Depth. The Three Curves on Each Graph are the Dependency of 
T| (m, n) on Time Lag nAr for m = 0 (no lag, dash-dot), m = 1 (10 km lag, solid), and 
m = 15 (150 km lag, dots). Data were Obtained from the MOODS Temperature Data 
Base for the Years 1929 to 1991. 
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Figure D.7 Spatial Variations of Autocorrelation Function t| (m, n ) for Type II at 
the Surface. The Three Curves on Each Graph are the Dependency of T| (w, n) on 
Spatial Lag (mAy) for n = 0 (no temporal lag, dash-dot), n = 1 (1 day lag, solid), 
and n = 15 (15 day lag, dots). Data were Obtained from the MOODS Temperature 
Data Base for the Years 1929 to 1991. 
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Figure D.8 Spatial Variations of Autocorrelation Function r| (m, n) for Type II at 
50% Water Depth. The Three Curves on Each Graph are the Dependency of 
r| (m, n) on Spatial Lag (mAy) for n = 0 (no temporal lag, dash-dot), n = 1 (1 day 
lag, solid), and n = 15 (15 day lag, dots). Data were Obtained from the MOODS 
Temperature Data Base for the Years 1929 to 1991. 
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Figure D.9 Spatial Variations of Autocorrelation Function t| (m, n) for Type 13 at 
80% Water Depth. The Three Curves on Each Graph are the Dependency of 
r| (m, n) on Spatial Lag (mAy) for n = 0 (no temporal lag, dash-dot), n = 1(1 day 
lag, solid), and n- 15 (15 day lag, dots). Data were Obtained from the MOODS 
Temperature Data Base for the Years 1929 to 1991. 
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Figure D.10 Spatial Variations of Autocorrelation Function rj (m, n) for Type in at 
the Surface. The Three Curves on Each Graph are the Dependency of rj (m, n) on 
Spatial Lag (mAy) for n = 0 (no temporal lag, dash-dot), n = 1 (1 day lag, solid), 
and n = 15 (15 day lag, dots). Data were Obtained from the MOODS Temperature 
Data Base for the Years 1929 to 1991. 
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Figure D.l 1 Spatial Variations of Autocorrelation Function r| (m, n) for Type IE at 
50% Water Depth. The Three Curves on Each Graph are the Dependency of 

r| (m, n) on Spatial Lag (mAy) for n = 0 (no temporal lag, dash-dot), n = 1 (1 day 
lag, solid), and n = 15 (15 day lag, dots). Data were Obtained from the MOODS 
Temperature Data Base for the Years 1929 to 1991. 
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Figure D.12 Spatial Variations of Autocorrelation Function t| (m, n) for Type IE at 
80% Water Depth. The Three Curves on Each Graph are the Dependency of 
t| (m, n) on Spatial Lag (mAy) for n = 0 (no temporal lag, dash-dot), n = 1 (1 day 
lag, solid), and n = 15 (15 day lag, dots). Data were Obtained from the MOODS 
Temperature Data Base for the Years 1929 to 1991. 
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